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Multigrid for propagators of staggered fermions in four-dimensional 
SU (2) gauge fields 
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Multigrid (MG) methods for the computation of propagators of staggered fermions in non-Abelian gauge fields 
are discussed. MG could work in principle in arbitrarily disordered systems. The practical variational MG methods 
tested so far with a "Laplacian choice" for the restriction operator are not competitive with the conjugate gradient 
algorithm on lattices up to 18 4 . Numerical results are presented for propagators in SU(2) gauge fields. 
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1. INTRODUCTION 

In Monte Carlo simulations of lattice gauge the- 
ories with fermions the most time-consuming part 
is the computation of the gauge field dependent 
fcrmion propagators. Great hopes to compute 
propagators without any critical slowing down 
(CSD) are attached to multigrid (MG) methods 
|I]-^0|. The MG methodology was reviewed by 
Brandt last year |Tl"| ]. 

We wish to solve an equation 



(-fl 2 + m 2 ) X = f 



(1) 



by MG, where is the gauge covariant staggered 
Dirac operator, and m is a small quark mass. 

The following MG notations will be used. The 
fundamental lattice of lattice spacing ao is de- 
noted by A . The first block lattice A 1 has lattice 
spacing a\ — L^aQ. Restriction and interpolation 
operators C and A, respectively, are given by ker- 
nels C(x,z) and A(z,x) with z e A , x e A 1 . If 
z is a site in block x, we write z £ x. Advantage 
is taken of the fact that we work in d = 4 di- 
mensions and that Lf, = 3 will be chosen, so that 
only a two-grid algorithm was implemented. The 
residual equation on the coarse grid was solved ex- 
actly by the conjugate gradient (CG) algorithm. 

2. BLOCKED STAGGERED FERMIONS 

We use a blocking procedure for staggered 
fermions which is consistent with the lattice sym- 

*Work supported by Deutsche Forschungsgemeinschaft. 



metries of free fermions 0. This forces us to 
choose Lb = 3. Even Lb are not allowed. 

Fig. [j] illustrates our choice of blocks. The 
different fermionic degrees of freedom are called 
"pseudoflavor" ||. Different pseudoflavors are 
distinguished by different symbols in Fig. [I]. 
Block centers x are encircled. The boundaries 
of seven blocks are marked. The averaging ker- 
nel C(x,z) is only nonvanishing if site z has the 
same pseudoflavor (symbol) as the block center x. 
Therefore the seemingly overlapping blocks have 
actually no sites in common. 

3. GROUND-STATE PROJECTION MG 

The averaging kernel C is chosen according 
to the ground-state projection (GSP) definition. 
The idea behind this definition is that the lowest 
mode of the fermion matrix, which is responsible 
for CSD in traditional algorithms, should be rep- 
resented as well as possible on the block lattice. 
GSP was introduced in Refs. @,||,|. An al- 
ternative MG approach to the inversion of the 
fcrmion matrix is the "parallel-transported MG" 
of Ben-Av et al. [00. 

In GSP the kernel C* of the adjoint of C fulfills 
a gauge covariant eigenvalue equation,^] 

(-A N , x C*)(z,x) = \ (x)C*(z,x) . (2) 

Xq(x) is the lowest (gauge invariant) eigenvalue 
of the positive (semi-)dcfinitc operator — An jX 

J C depends on the gauge field although this is not indi- 
cated explicitly. 
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Figure 1. Blocking of a two-dimensional stag- 
gered lattice consistent with the symmetries. 



which will be specified below. C*(z,x) = 
C(x, z)\ where t denotes Hermitean conjugation, 
as usual. The solution of Eq. (||) is made unique 
by imposing the normalization condition 

CC* = 1 (3) 

and the covariance condition 

C(x, x ) = r(x) I (4) 

with r(x) > 0. (In case of gauge groups Q differ- 
ent from £7(1) or SU(2) the r.h.s. of is replaced 
by a positive Hermitean matrix.) The covari- 
ance condition ([|) ensures that the kernel C{x, z) 
transforms covariantly under gauge transforma- 
tions. In the limiting case of vanishing gauge cou- 
pling, GSP yields (gauge-transformed) piecewise 
constant kernels on blocks. 

In a gauge theory with N c colors, C(x, z) is an 
N c x N c matrix which will not be an element of Q 
in general. For Q = SU (2) it is easy to see that 
1-column vector solutions of (^J) may be combined 
into 2-column matrix solutions [pj . 

3.1. Bosons 

In case of bosons Ajv jX is the covariant Lapla- 
cian with Neumann boundary conditions on 



boundaries of blocks x. Ajy _ x is defined by 

{A NjX C*)(z,x) = 
J2 [U(z,z')C*(z',x)-C*(z,x)} (5) 

z n.n. z 

z' £x 

for z e x. U(z,z') is the gauge field on the 
link (z, z'). As a supplement we enforce Dirichlet 
boundary conditions in such a way that C(x, z) 
unless z G x. 

3.2. Staggered fermions 

For fermions A in (||) was chosen as the gauge 
covariant fermionic "2-link Laplacian" which is 
defined through -jp 2 = -A + a^F^ §). With 
this "Laplacian choice" C(x, z) is only nonvan- 
ishing when z and x carry the same pseudoflavor. 
Other proposals were made in Ref. |(|, but have 
not been implemented yet. 

There exists an efficient algorithm for the solu- 
tion of Eq. (|) 01. Therefore it is not necessary 
to sacrifice the gauge covariance of GSP by wast- 
ing computer time for gauge fixing. 

3.3. Idealized MG algorithm 

Given the averaging kernel C there exists an 
associated idealized interpolation kernel A. By 
means of this A a (nearly) critical system on A 
is mapped onto a noncritical system on the MG. 
The construction of the optimal A was described 
in Mack's Cargese lectures [^5|. Its origin are 
works on constructive quantum field theory. 

The gauge covariant optimal A is the solution 
of the equation 

([D + kC*C]A)(z,x) = kC*(z,x) (6) 

for large k. D = —A + m 2 for bosons, and D = 
—p 2 + m 2 for fermions. 

An averaging kernel C is said to be "good" 
(or to define a "good" block spin) if the associ- 
ated A(z, x) decays exponentially in \z — x\ . GSP 
defines a good C in arbitrarily disordered gauge 
fields, both for bosons and for staggered fermions. 
Fig. H visualizes an example. A decays exponen- 
tially over distance L^a^ (nearly) as fast in the 
presence of gauge fields than in their absence. 

By means of the idealized A it was shown in 
Ref. (8| that MG computations of bosonic prop- 
agators without CSD are possible in arbitrarily 
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Figure 2. An optimal fermionic interpolation ker- 
nel A(z,0) in a quenched SU(2) gauge field at 
j3 = 2.7 on an 18 4 lattice. Shown is a two- 
dimensional cut through the block center x = 0; 
Z3 and 24 are fixed. The vertical axis gives the 
trace norm of A{z, x). 



disordered gauge fields (including the case j3 = 0). 
These computations showed for the first time that 
the MG method could work in principle in arbi- 
trarily strong disorder. Unfortunately, the ideal- 
ized MG algorithm is not practical for produc- 
tion runs because of computational complexity 
and storage space requirements, but it was im- 
portant to answer questions of principle. 

3.4. Variational MG 

A practical MG method is variational coarsen- 
ing. There A = C* , and the coarse grid operator 
is CDC* . This method reduces to piecewise con- 
stant interpolation in the absence of gauge fields, 
where it is known that CSD is eliminated. 

However, in nontrivial gauge fields variational 
MG with the Laplacian choice of C is not com- 
petitive with CG on lattices up to 18 4 ; see Fig. |[ 
The non-optimized MG program requires a factor 
of 4.5 more arithmetic operations than CG. 

A further numerical finding is the following. 
Relaxation times r of conventional one-grid relax- 
ation and of variational MG with the Laplacian 



Figure 3. Computation of propagators of stag- 
gered fermions (— Tp + to 2 ) -1 on an 18 4 lattice 
in a quenched SU{2) gauge field at (5 — 2.7. 
r denotes the residual. The numbers refer to the 
following algorithms: 1/2: variational MG SOR 
(lexicographic ordering) with to 2 = 0.01/0.0, u) = 
1.88/1.92; 3/4: CG with to 2 = 0.01/0.0. 



choice of C follow a scaling relationg 

COTlst. , 1 a? 2 2 

t = — — 5- with Am = to — to , (7) 

where to 2 ,, is the lowest eigenvalue of — p 2 , and 
const, is independent of the lattice size. m 2 r is 
small and is usually neglected, but on relatively 
small lattices this neglect is not justified. From 
([?]) we conclude that variational MG with the 
Laplacian choice of C will not be able to elim- 
inate CSD on large lattices. 

4. UPDATING ON AN MG LAYER 
CONSISTING OF A SINGLE SITE 

Because of (Q), (MG) relaxation continues to 
have CSD even on a lattice of only 2 d sites. 
Therefore it seems necessary to update on a "last 
site" in order to eliminate the appearance of CSD. 

Given an approximate propagator x^ n ' , this up- 
dating (in the unigrid point of view) amounts 
to globally rescaling x by a gauge covariant 
N c x N c matrix f2 (in case of bosons or Wilson 

2 In case of one-grid relaxation for bosons Eq. {^) is known 
analytically. 
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fermions) fllRj: 

X (n) (z)^X (n) (z)fl . (8) 

The matrix Q is chosen such that the energy func- 
tional K[x] = \ <x, D X > - < x, / > of the 
rescaled approximation X ^ Q is minimized. 

In case of staggered fermions we have to con- 
sider that there are 2 d different pseudofiavors, and 
we replace ft in (0) by Q(H(z)) where H(z) de- 
notes the pseudoflavor of z. The 2 d N c x N c matri- 
ces £l(H) are again chosen automatically by the 
algorithm in such a way that K[x^&\ gets min- 
imized. 

We note that Q, and the Q(H)'s practically 
equal 11 as soon as errors decay exponentially. 

4.1. Results for bosons 

In a fixed volume the asymptotic fall-off prop- 
erties do not depend on Am 2 when (^) is in- 
cluded in algorithms. When the lattice size is 
increased, the asymptotic fall-off rate does not 
change, there is only a mild volume effect with 
respect to how long it takes until errors decay 
exponentially. This volume effect is much milder 
than the one in (unpreconditioned) CG, and vari- 
ational MG plus (@) begins to become superior in 
computer time for lattices >18 4 . 

4.2. Results for staggered fermions 

In a fixed volume the asymptotic 1 / Am 2 di- 
vergence is eliminated when the analog of (JsJ) is 
included in algorithms. However, it takes a rela- 
tively large number of iterations until fall-offs are 
exponential. Also, in contrast to the bosonic case, 
asymptotic decay rates increase at fixed Am 2 
when the lattice size is increased. Therefore CSD 
seems not to be eliminated, and it is questionable 
whether the method pays in practice on lattices 
of realizable sizes. For details see Refs. |p|,^6|. 

5. CONCLUSIONS 

The variational MG methods tested so far for 
staggered fermions are not competitive with CG. 
We expect better results when blocks are used 
which overlap in a nontrivial way. The question, 
whether it is sufficient to average only over sites 
with the same pseudoflavor in nontrivial gauge 



fields, also has to be investigated. The Diracian 
proposal for C of Ref. § should be tested. 
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